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ABSTRACT 


A theory is presented to explain the observed variations of the magnetic 
field and plasma in the vicinity of the moon. Under the guiding-center 
approximation, solutions for the plasma flow near the moon are obtained from 
the kinetic equation. The creation of a plasma cavity in the core region of 
the lunar optical shadow disturbs the interplanetary magnetic field. Maxwell's 
equations are used to study perturbations of the magnetic field in the lunar wake. 

The acceleration drift current, which was omitted from the earlier work, is 
included in the present theory in the calculation of the total electric current 
in the lunar wake. Numerical solutions of Maxwell's equations are obtained. 

When the interplanetary magnetic field lines penetrate into the lunar body, due 
to sudden change of magnetic permeability the magnetic field is disturbed at the 
lunar limbs. Propagations of this disturbance with magnetoacoustic speed form a 
Mach cone downstream, which is sometimes observed as the exterior increase^ of 
field magnitude in the lunar penumbra. Perturbations of the magnetic field are 
restricted to the region inside the Mach cone; the region outside remains undisturbed. 
The numerical results agree extremely well with experimental data from the 
Explorer 35 spacecraft. 


I. INTRODUCTION 


Measurements^ of the interplanetary magnetic field and plasma in the 
vicinity of the moon have been made from lunar orbit on the Explorer 35 spacecraft. 

The purpose of this paper is to present a theory which can explain the observed 
variations of the field and plasma in the lunar wake. 

When the solar wind interacts with the moon, no shocks are observed in 
the vicinity of the moon. Figure 1 shows a simultaneous measurement of the 
interplanetary field and plasma on Explorer 35 when the moon is outside the earth's 
bow shock. The major effect of the moon on the solar wind plasma is the creation 
of a plasma cavity in the umbra I region of the lunar shadow. In this cavity the 
magnitude of the magnetic field increases, i.e., the field is observed to be 
stronger than the undisturbed interplanetary field. On either side of the umbral 
increase, the field decreases, ?.e 0 , the field becomes weaker than the undisturbed 
condition. These penumbra! decreases occur at the location where the plasma 
density is about half of the undisturbed plasma density, and they are often bounded 
on the exterior by additional small increases in the field magnitude. A positive 
correlation between the increase of the plasma flux, the plasma density and the 
penumbra I increase of field magnitude has also been observed. 

From the kinetic theory point of view, the plasma flow in the vicinity of the 
moon can be analyzed under a guiding-center approximation. The zeroth-order 
solution for the plasma flow obtained by Whang^ is in good agreement with the measured 
plasma density, plasma flux, and plasma flow direction . 
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The first order guiding-center theory is used to describe the plasma and 

O 

the electric current in the lunar wake. Whang 0 has previously studied the 
perturbation of the magnetic field in the vicinity of the moon by assuming that 
the total electric current in the lunar wake is composed of the gradient drift 
current, the curvature drift current and the magnetization current. This present 
theory produces more realistic results than the earlier theory in three aspects: (1) 

The present theory infers that perturbations of the magnetic field in the lunar 
wake are confined to a region bounded by a Mach cone (a standing magnetoacoustic 
wave), outside of which the interplanetary field is undisturbed. (2) The penumbral 
decreases of the field magnitude are of the same amplitude as the observed decreases. 
(3) The umbra I increase of the field magnitude is of the same spatial extent as the 
observed increase. The basic difference of the present theory from Whang's previous 
one is the inclusion of the acceleration drift current in the calculation of the total 
electric current in the lunar wake. This current turns out to be at least as large 
as the other terms. 

The solar wind flows approximately radially outward from the sun at 1 AU 
at a supersonic velocity. The solar magnetic lines of force are carried outward by the 
solar wind and are twisted into the form of Archimedean spirals, thus the interplanetary 
magnetic field is in general not aligned with the direction of the solar wind velocity. 
Thermal motion of the solar wind- ions is anisotropic; the ratio of T n and T x varies 
between 1 and 4. In the umbral region immediately behind the moon, only those 
ion particles with high parallel thermal energy are present, while the perpendicular 
thermal energy is directly proportional to field intensity. Thus, the plasma thermal 
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anisotropy, Tjj /T ± , increases in the core of the lunar wake. The propagation 
speed of magnetoacoustic waves in an anisotropic plasma is again highly anisotropic 
when the thermal energy of the anisotropic solar wind plasma is of the same order 
of magnitude as the energy of the interplanetary magnetic field. All these reul 
physical conditions of the solar wind are taken into account in the present model. 

The solar wind interaction with the moon has also been studied by other 
authors^'^"^ using the continuum theory approach. In these models, the thermal 
anisotropy of the solar wind plasma and the anisotropic acoustic waves are not 
considered, and a simple polytropic relation for pressure and density is used. 

Some have only considered the case where the magnetic field upstream of the moon 
is aligned with the direction of the solar wind velocity. 

Two dimensional steady-state solutions of Maxwell's equations are studied. 

Thus the present results actually correspond to the interaction of the solar wind 

with a cylindrical moon. Under the two-dimensional cylindrical-moon approximation , 

the solution will become less accurate at large distances behind the moon. 

The first order solution of the kinetic equation is also studied in the last 


part (section VII) of this paper. 
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II. KINETIC EQUATION 

The Larmor radius of the ion particles in the solar wind is only a small 
fraction of the radius of the moon, therefore, the lowest-order solution for the 
ion flow in the vicinity of the moon can be calculated^ by approximating each 
ion particle by a guiding-center particle. It is convenient to use the guiding- 
center distribution function^ which is defined as a function of x , t, v and p 

P v, (l) 


where v is the component of the particles velocity parallel to B, and p the magnetic 

13 

moment. From the first order theory of guiding-center motion , 






u x 


( 2 ) 


and 


— V- = O 

a-t 


( 3 ) 


where d/d* = 5M + ( Ux + & v) • V , $>=&/!}, Ux = C ExB/fi 2 

and 7< s : ' e, - ^7 # 

The collisionless Boltzmann equation can be expanded following the Chew- 
Go Idberger-Low theory^' The lowest-order equation can be reduced to 


«©(F/B) * O 9 (4) 

where the operator X) denotes the time derivative following the trajectory of a guiding- 
center in phase space. 



Equation (4) states that F/B is conserved following the trajectory of a guiding-center 
in phase space. 

Solutions of the kinetic equation (4) are studied subjected to the boundary 
conditions: (?) the moon absorbs all particles that hit its surface, and (ii), the 
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distribution function is Maxwellian in the undisturbed flow upstream. Once F is 


known. 

the ion density. 

parallel 

velocity, and parallel 

pressure can be calculated 

from 






" =1 


dir 

; 

(6) 


U.- ■ 


f dpt dir 

(7) 

and 

?»=( 

f rn. ( v 

2 

- Ufi) P djudu^ 

(8) 


The average magnetic moment is related to the perpendicular pressure by 



Consider that in the undisturbed upstream flow, the solar wind moves with a 
constant velocity U Q with respect to the center of the moon, and carries with it a 
steady uniform magnetic field B Q . We choose a coordinate system (Fig. 2) with its 

' mmmi 

origin located at the center of the moon, the x axis parallel to U 0 , and the z axis 
along the direction of the electric field in the upstream E o — - U© * § c / c . 
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The scale of the coordinate system is normalized by using the radius of the moon 
as a unit length. 

The zeroth order equations are those in which the guiding-center velocity 
of each particle is assumed to be constant in both magnitude and direction. In 
this order the equation may be solved analytically. The earlier theory' 7 presented 
the solution for the special case in which the magnetic field was constant. Here 
an analytical solution is presented in which the magnitude of the field is allowed 
to vary. 

As in the earlier work the space near the moon is divided into four regions 
as shown in Fig. 2. Let <}> 0 denote the direction angle (the angle between -U 0 
and B 0 ) in the upstream region. The analytical zeroth-order solution gives 

= P„/P„o = B/B 0 , 


and 


U„ = - U, «* 4> a 


in region 0. In the other regions the solutions are functions of B/B 0 , the two solar 

</j« 

wind parameters 4> o and S- [S=U 0 /< zkT tlo /fn.) J^and the two position parameters 


X = arcsin [(l - z 2 )' / 3 (x 2 + y*f l/2 ] , 


and 


The solutions are : 


A = arc-fan ( y/ x) t 


7T, = 2 % [ « rf c ( Y, S ) + erfc ( V” 2 S^] 


( 10 ) 
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, «*PC-y/S^e«pC-vW (11) 

^0 ° 2TT 1 ' 1 S nBo/”'B 

and 

•— -z=(n/n 0 )[ 1 - S 2 C co * <& + U„/U 0 f] + 

MIO 

+ (Tr? /i C8/8 0 )S[v l mMpC-YSSV+ytWC-tfS 9 )]' (12) 

where 

Y, ^ sin (\ + \) / i i n C c £©”X’"^) in regions 1,3 

— 00 in region 2 

V 2 = tin (\-X)A; n (4>.*X-\) in regions 2,3 

= 00 in region 1 , 

Equations (10) - (12) are slightly different from the solution given in reference 
7, because the frozen-in effects of the variation of field magnitude are included in 
the present form. With the inclusion of the factor B/Bo, the present solution is capable 
of explaining the positive correlation between plasma density and magnetic field 
strength in the exterior penumbra. 

A solution to the first order equations in which the deflection of the guiding- 
center trajectory in the non-uniform magnetic and electric fields is taken into account , 
»S‘ obtained numerically in section VII. 
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III* MAXWELL'S EQUATIONS 

Maxwell's equations are used to describe the variation of the magnetic 
field. The total electric current is composed of the magnetization current^ the 
gradient drift current, the curvature drift current and the acceleration drift current, 
Let D/Dt denote the derivative 2/b-fc + U • V , U the fluid velocity and p the 
fluid density. The acceleration drift current density is 

DU 


J * 4 e,* 

B *** ' Dt 


The acceleration drift current is sometimes known as the polarization current^, and 
was not included in the early theoretical treatment^/S of the solar wind interaction 
with the moon. This term is at least as important as the other terms of the current 
density. Maxwell's equations can now be written as 

V B - O (13) 


and 

( 1 + -^pr) Vx 8 * Pj. vin n * P„ X + f (14) 

Here n is the number density of the plasma, P x and P„ respectively the perpendicular 
and the parallel plasma pressure. 

In the wake region, the electric field is only slightly disturbed from the 
undisturbed electric field E 0 . The perpendicular component of the fluid velocity can 
be calculated from 
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Uj. = c E„ * B / B z (15) 

The perpendicular pressure is directly proportional to the product nB. Once n,' U„ 
and P (| are known functions of x and B, Equations (13) and (14) can be solved 
for the variation of B in the lunar wake. 

Two dimensional steady-state solutions of Maxwell's equations are studied 
in this paper. 


IV. THE EQUATIONS OF CHARACTERISTICS 


In a collisipnless plasma, the magnetoacoustic wave^ propagates along 

1/2 

the field direction with speed ci=s[(P x -P (l +-B/^7r)/f] , and perpendicular to the 

^ _ j /2 

field direction with speed b=[(2Pj r +B >*V it)/^ J . The solar plasma flow is 
supersonic in the vicinity of the moon except in the umbral region of the lunar 
wake, where the magnetoacoustic speed becomes very large because the plasma 
density and pressure decreases. In the supersonic region outside the lunar umbra 
the method of characteristics is used to calculate the variation of the magnetic 
field. 

let B and 4> denote respectively the magnitude and the direction angle 
of the vector JB, and ®2 = e z x e^ , then one can write 

v*B = (- e* • < 7 b + B f i * ^ 4 > ) e x 

6, k >c s <7 x e, 

= C$i • ^40 € z 

and r U *(7U A e» + U„ 7e f )x §, 

= U'C-U^ vln B + U||^4>)«x . 

Making use of these relations, one can write Maxwell's equations {2 (13), (14)] 

£» • V in B + £* * = ° 

and A , • V in B + A* • V <#> = H 

where 
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A 2 = o 2 e, - U„ U / 

and M - C Rl /p) £2 * V in (n/B)' 

We now consider the linear combination of the two equations (16) and (17) 
obtained by multiplying equation (16) by a factor t and add the two equations. 


L • V In B 4 - M * V<f> ~ H 


where 


L = C u± u„ + -c ) e, + C 2 - b 2 ) e 


and 


M = ( a 1 - u*) e, + C •* ~ U* U„ ) e* 


(18) 


The magnitude and the direction of vectors L and M^may be adjusted by the variation ’ 
of the factor t. Now, suppose for some choice of t, the vectors L and iM are parallel 
and tangent to a curve C; then only the derivatives along the C-direction are contained 
in equation (18). In this case the curve C is known as a characteristic curve. 

L and M are parallel vectors when 

T 2 = a 2 Uj+ b 2 ( V - « 2 *> 2 

or 

t + = -T- - («»*«/ + b 2 u „ 2 - b* )' * 

Substituting these two roots into (18), one can replace the Maxwell's equations by 

a pair of characteristic equations 


L + • V bn B 4 M + * ^ = H 

L • v in 8 -4 M = H 
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where 

L± = S' + At 

M *— > 

M. = *r A . e 2 + A 2 

L. and hA are parallel vectors, thus each of the above equations contains the 
derivative of B and <f> along the direction of a characteristic curve. 

The characteristic equations (J9) are used to calculate the variation of 
magnetic field in the supersonic region (where > 0) of the lunar wake. In the 
subsonic and transonic region a different method, which is developed in the next 
section, is used in calculations. 

In equation (19), if all three terms in one of the two equations are zero, 
its solution may be called the simple wave solution in frozen-in plasmas. 
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V. THE UMBRAL ANALYSIS 

The method of characteristics derived in section IV can be used to calculate 
the variations of magnetic field in the supersonic region of the lunar wake. In the 
subsonic and transonic regions of the lunar wake, approximate solutions of Maxwell's 
equations are calculated under the assumptions (i) the magnetic field, J3, changes 
rapidly in the y-direction, but very slowly along the x-direction, and (?i) B x and 
By are of the same order of magnitude. 

The divergence of _B equation requires that 2> B*/dx = ~2>B y /by # Thus 
the order of magnitude for the derivatives of can be measured by 


IBx . 2> B x , B y . lb By 

’ Sy ‘ ax 


1 : e : € : € 


Dropping out the second order terms, Maxwell's equations (13), (14) may be 
written as 


7 -7 + ? sB y 


( 20 ) 


where 


Z,= ct 2 B/ + b* B. 2 - U 2 B 


2 2 * 


Z 2 = - ( B 3 Pi /p) e, x \7Ln C n /B) 
and £ 3 = 2 ( a 2 - b 2 ) B x B y -2 B 2 U x U y . 

For given x, the method of characteristics is used to determine the solution 
for B in the supersonic region outside the lunar umbra. Thus, conditions at the lower 
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and the upper boundary of the umbra are determined from the characteristics 
solution. 


B x B x ] 

at 

y = '/] 

(21) 

Bx ~ B X 2 

at 

y = y 2 

(22) 

By = B yl 

at 

y = y\ 

(23) 

B y = By 2 

at 

y = y 2 . 

(24) 


The divergence of B equation also requires that, when x is fixed 

y i 

+ 

y, 

Subjected to the above five conditions, Eq. (20) can be integrated numerically for 
the solution of B in the interval y^ <£ y ^ y 2 » 

In Eq„ (20), the density gradient in Z 2 is the main driving force for the 
variation of B x along the y-direction in the umbral region a Integration of Eq. (20) 
is carried out by assuming a series form for the first order term d B y /<$y » Let 

By = c 0 + c, y + c, y* + c, y 3 ^ ( 26 ) 

where the four constants c Q , C], C 2 , C 3 are adjusted to satisfy the conditions (21)-(25) 0 

Making use of the condition (23) and (24), C 2 and c^ can be determined as functions 

of c 0 , C|. As the first order term is a known function of c Q , cj and y, 

numerical integration of Eq. (20) can be carried out starting at y = y € with B = B vl . 

X x x 1 

The solution for B x will be in the form of 


)( ^5 4> a )c/y = o # (25) 



-Oo 
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B„ = B* ( c«, <=i, y) (27) 

Here we still have two adjustable constants c Q and c^ , Their values can be determined 
when both Eq. (22) and Eq 0 (25) are satisfied . The method described above indicates 
that by assuming a cubic polynomial solution for the first order term, one can calculate 
the solution for B x satisfying all required conditions,. 
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v| . NUMERICAL SOLUTION 

Given the speed ratio S and direction angle <p 0 , the zeroth order solutions 
of n/n 0 , Ujj /U Q and P„ /P, lo for ion particles can be calculated from Equations 
(10) - (12). These analytical solutions are used to approximate the plasma conditions 
in the lunar wake. Then based on the mathematical method developed in the two 
preceding sections, solutions for the perturbation of magnetic field in the lunar 
wake can be carried out numerically. The results quantitatively reproduce the 
major features of the observed magnetic field perturbations in the lunar wake. 

The perturbations of the magnetic field in the lunar wake depend upon four 

solar wind parameters: the plasma anisotropy T lt0 /T xo , the speed ratio S, the 

direction angle and the parameter pafctrPxo/B 2 . G is defined as the ratio 

of the perpendicular thermal pressure to the field pressure (the plasma thermal 

pressure includes both the ion and the electron thermal pressure). For typical 

interplanetary conditions at the orbit of the earth, the ratio of T ,j and T x varies 

between 1 and 4, the speed ratio S is of the order of 7^10, the magnetic field 

has a direction angle ~ 135°or 315° and the G value is of the order of unity. 

Experimental evidence has shown that the field anomalies are strongly affected 
18 

by the G value 

A numerical solution for some typical solar wind parameters is plotted in 
Figures 3 and 4. 

Figure 3 shows the variation of the field magnitude in the lunar wake region . 
The exterior increase of the field magnitude actually represents a standing magneto- 
acoustic wave attached to the lunar limbs; it resembles the Mach cone in aerodynamics. 



Perturbations of the magnetic field are restricted to the region bounded by the 
Mach cone; the region outside remains unperturbed. 

Figure 4 shows that the width of the umbral increase is always slightly 
less than one lunar diameter. This feature agrees with observed data, and was 
not predicted by previous theories. The penumbral decrease of field magnitude 

8 

obtained here again agrees with measurements far better than Whang's early result . 

The numerical solution is carried out by assuming that the field is slightly 
disturbed at the initial line x = 0 near the lunar limbs. These disturbances propagate 
downstream to form the standing magnetoacoustic wave, which is observed as 
penumbral increases, if the field is completely undisturbed at the initial line, 
then numerical solutions cannot produce the observed penumbral increases. This 
suggests that when the magnetic field lines penetrate into the lunar body, due to 
sudden change of magnetic permeability across the sunny side of the lunar surface, 
the magnetic field is disturbed near the lunar limbs. These disturbances propagate 
in the supersonic flow to the downstream. If the field magnitude increases at the 

t 

limbs, the propagation of this disturbance would be observed as a penumbral increase 
of the field magnitude in the lunar wake. 
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VII. FIRST ORDER SOLUTION OF KINETIC EQUATION 

The zeroth order solutions of the kinetic equation discussed in section II 
are obtained under the assumption that the guiding-center particles are not 
accelerated; their trajectories are approximated by straight lines. However 
since the electric and the magnetic fields are perturbed in the lunar wake, the 
velocity and the trajectory of a guiding-center will also be perturbed. The effect 
of the induced electric field on the guiding-center trajectories has been discussed 
by Alfonso-faus and Kellogg^. In this section the combined effect of the induced 
electric field and the perturbed magnetic field are included in calculating the 
first order solution of the guiding-center kinetic equation. 

We consider that in the perturbed region the electron density is equal to 
the ion density and the electrons always assume an equilibrium distribution in the 
potential field. The potential of induced electric field depends on the ratio n/n 0 
according to the logarithmic law 

^ = ( feT t /e) L(n/n 0 ) (28) 

From this equation, one can calculate that the magnitude of the induced electric 
field is very small compared with the upstream electric field E 0 . 

The guiding-center kinetic equation states that F/B is constant on the path 
of a guiding-center in phase space 


8 ( F/ B) = ° 
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One may write the operator JQ in the following form 

A = R.& + R z + 

which indicates that the guiding-center trajectory is described by 

dx _ dy d v 

R, R 2 R* , 


(29) 


(30) 


The expressions for Rj's are of the following forms; 


R,= Uj. ^ “ v cos , 

R 2 = - c«5 4* - IT sm d> ; 

and ^=5T= - e , + 

* 

The procedure of calculating the first order solution for the distribution 
function F(v,p) at any point in the perturbed region is the following. If B and ^ 
are known, then one can integrate Eq. (30) to obtain the trajectory of a guiding 
center and the variation of v a long its trajectory. The magnetic moment p is an 
adiabatic invariant as the guiding center moves. In the undisturbed region upstream, 
the guiding-center trajectories are straight lines, and the distribution function is 
Maxwellian, F^. If the trajectory intercepts the lunar surface, the distribution 
function F(v,p) is zero. Otherwise, the trajectory will end at the undisturbed 
region with its parallel velocity changed to v 0 , the distribution function can now be 
calculated from 


F(v, ^0 = ( B/ B») F.O®, fO. 


(31) 
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Following this piocedure, one can numerically calculate the- guiding-center 
distribution function at any point in the perturbed region as functions of vand p. 
Taking moments of this distribution function one can calculate the first order solution 
for n, Uj| and P„ at the point of interest. 

One further simplifying assumption is used to carry out some first order 
solutions in this paper, that is, all particles are assumed to possess the same 
magnetic moment 

^ “ PiO / "o , 

Under this assumption, the guiding-center plasma has thermal motion only along 
the field line; its distribution function depends on one velocity variable V . The 
plasma behaves as a one-dimensional gas. 

Consider that (i) the electron temperature Te = 2T no , (ii) the potential 
is calculated from the zeroth-order n, and (iii) the magnetic field JB is calculated 
from the method outlined in the previous sections. One can calculate the first- 
order solution for F. Some first order solutions are plotted in Figure 5. The first 
order solution shows a slightly decrease of the plasma density gradient across the 
boundary. The analytical solution for the zeroth order density is practically good 
enough for the purposes of understanding the plasma flow in the lunar wake. 

Figure 5 also shows that (i) the penumbral decreases of field magnitude occur 
at the location where the plasma density is about half of the upstream density, and 
(ii) a positive correlation exists between the increases of plasma density and magnetic 
field strength in the penumbra. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 
Figure 3 


Figure 4 
Figure 5 


The simultaneous measurements of field and plasma obtained on 
August 5, 1967 from lunar orbit on the Explorer 35 spacecraft. 

The trajectory of the spacecraft is shown projected on the 
ecliptic plane and positionally correlated with the data through 
UT annotation. The x axis is parallel to the sun-moon line. 

Four characteristic regions of perturbed plasma flow in the vicinity 
of the moon. 

Distributions of the magnitude of magnetic field in the lunar wake. 
Perturbations of the field are restricted to the region bounded by 
the Mach cone. The dashed lines indicate the direction of plasma 
velocity. 

Perturbations of the field magnitude in the lunar wake calculated 
for T llo /Ti e =M.5, S=io, <*> o = \35°, £=1. 

The zeroth order and the first order solutions of the ion density. The 
first order solutions are calculated for Te = 2T uo . 
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